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Abstract 

We investigated the conductance and shot noise properties of qnasiparticle-transport through 
a superconducting barrier in graphene. Based on the Blonder, Tinkham, and Klapwijk (BTK) 
formnlation, the theory to investigate the transport properties in the superconductive graphene 
is developed. In comparison, we considered the two cases that are the transport in the presence 
and absence of the specular Andreev reflection. It is shown that the conductance and shot noise 
exhibit essentially different features in the two cases. It is fonnd that the shot noise is suppressed 
as a result of more tunneling channels contributing to the transport when the superconducting 
gate is applied. The dependence of the shot noise behavior on the potential strength and the width 
of the superconducting barrier differs in the two cases. In the presence of the specular Andreev 
reflection, the shot noise spectrum is more sensitive to the potential strength and the width of the 
superconducting barrier. In both of the two cases, total transmission occurs at a certain parameter 
setting, which contributes greatly to the conductance and suppresses the shot noise at the same 
time. 
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I. INTRODUCTION 


In the past decade, there has been a great deal of interest in stndying the physical 
properties of graphene both theoretically and experinientallyi“-. In graphene, the low- 
energy excitations are massless and chiral Dirac fermions with linear energy dispersion near 
the Dirac point. Because of the unique energy dispersion in graphene, there are some 
special phenomena such as the unconventional quantum Hall effect, Klein paradox, eto^^-*^. 
Particularly relevant to the present work, the presence of the Dirac Fermions in the graphene- 
based superconductor junctions results in the specular Andreev reflect ioni*^. 

The Andreev reflection^ describes the tunneling phenomenon of electron excitation con¬ 
verting into hole excitation by the superconducting pair potential. The interface between a 
metal and a superconductor can reflect a negatively charged electron incident from the metal 
into a positively charged hole, while the missing charge enters the superconductor forming a 
Cooper pair. Usually the hole is reflected back along the path of the incident electron in the 
conventional materials, which is called “retro-Andreev reflection”. However, the Andreev 
reflection in undoped graphene is specular—the so-called ’’specular Andreev reflection”, in 
which the reflected angle is invertedi*^*^*^. Because the graphene needs to be described 
by the Dirac-like equation rather than the usual Schrodinger equation, the superconduc¬ 
tive graphene needs to be described by the Dirac-Bogoliubov-de Gennes equation rather 
than the Bogoliubov-de Gennes equation for usual superconductorsi‘^'ii'l^'l^>^. However, 
graphene is not a natural superconductor. Recent research has shown that superconductiv¬ 
ity can be induced in a graphene layer in the presence of a superconducting electrode by 
means of the proximity effectf^'^'^. So far, the physical properties of the graphene-based su¬ 
perconductor junctions have been extensively studied and many important results have been 
obtained. Most of them focus on the transport properties of the graphene/superconductive 
graphene interface and only the conductance is consideredi'^'^^^— , which prompts us to 
investigate the transport properties through a superconducting barrier in graphene and es¬ 
pecially focus on the shot noise spectrum. In this work, based on the Dirac-Bogoliubov-de 
Gennes equation and the Blonder, Tinkham, and Klapwijk (BTK) formulation, we derive 
the shot noise formula and provide numerical results of the transport properties in the 
graphene/superconductive graphene/graphene (G/SG/G) heterostructure. Dependence of 
the shot noise spectrum on the structure parameters and its physical mechanisms are dis- 
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cussed. 


II. THEORY AND MODEL 


We consider transmission through the G/SG/G heterostructure occupying the x-y plane, 
the schematic of which is shown in Fig. 1. The growth direction of the graphene is taken 
to be the a:-axis. The left G region extends from x = —cxo to a: = 0 and the right G region 
extends from x = a to x = +oo. The superconductive region occupies the 0 < x < a 
region. The electron and hole excitations are described by the Dirac-Bogoliubov-de Gennes 
equationi'^'ii^— : 


H^-Ef a (r) 
A*(r) Ep-H^ 




( 1 ) 


Here, = {'ipAa,'ipBa,'ipAay ~'^Ba) four-component wave functions for the electron 

and hole spinors. The index a denotes K or K' for the electrons or holes near the Dirac K 
and K' points, a takes values K' [K) for a = K {K'). Ep denotes the Fermi energy. A and 
B denote the two inequivalent sites in the hexagonal lattice of graphene. The Hamiltonian 
Ha is given by 


Ha = -ihvp [a^d^ + sgn (a) aydy] + U (r), 


( 2 ) 


where vp denotes the Fermi velocity of the quasiparticles in graphene and sgn (a) takes the 
values of ±1 for a = iF (K'). 

For 0 < X < a, the superconducting electrode on the top of the graphene layer induces a 
nonzero pair potential A (r) via the proximity effect. We model the pair potential as 


A(r) = { 


0 , 

Aoe'^ 


others, 

0 < X < a. 


( 3 ) 


where Aq and <p are the amplitude and the phase of the induced superconductive order 
parameter. The electrostatic potential U (r) in the G and SG regions can be tuned indepen¬ 
dently by a gate voltage or by doping. We take 


G(r) = { 


0 , 

-Go, 


others, 

0 < X < a. 


( 4 ) 


Eq. ([!]) can be solved straightforwardly to yield the wave functions T in the G and 
SG regions, respectively. In the G region, for the electrons and holes traveling in the ±x 
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directions with a transverse momentum ky = q and energy e, the wave functions are given 
by 

^ exp {^ia/2) ^ 


e± _ exp(iqy^ikx) 


'^N = 


y/cosa 


/ 


h± exp(iqy±ik'x) 

Vcos a' 


^ N — 


±exp {±ia/2) 

0 
0 
0 
0 

exp (=Fia'/2) 
y Texp{±ia'/2) j 


( 5 ) 


with 


»■"«' = 51:. 

^ _ e+Ep 


( 6 ) 


T-^ COS a, 
avp ’ 


k' = ^ 7 -^ cos a'. 

hvp 

a is the incident angle of the electron and a' is the reflection angle of the hole. Note that 
for the Andreev process to take place, the maximum incident angle of the electron is given 
by 


ttc = arcsm --— . 17) 

\£ + Ef J ^ ^ 

In the SG region, the quasiparticles are mixtures of the electrons and holes. The wave 
functions of these quasiparticles moving along the ±a;-direction with the transverse momen¬ 
tum q and energy e has the form of 


/ 


^e± _ ^iqy±i{ko-iK) 


\]>^± — ^iqyTi{ko+iK.) 


exp {—i/3) 

±exp {—i/3 ± iy) 
exp {—i(/>) 

±exp (iiy — i(/>) J 
( exp {i/3) ^ 

^exp {i/3 =F iy) 
exp {—i(/>) 

^exp {Til — i(/>) J 


( 8 ) 
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The parameters /?, 7 , ko, and k are defined by 


arccos {e/ Aq) , e < Aq, 

L 

—iarcosh {e/ Aq) , £ > Aq, 

7 = arcsin [hvpq/ {Uq + Ep )], 
kQ = sJ{Ep + UQf/{hvpf-q\ 

in^ 

{hvpYko 


( 9 ) 


Taking into account both the Andreev and normal refiection processes, the wave functions 
in the left G, SG, and right G regions can be written as 


^2 = A%+ + B%- + GT^+ + (10) 

^3 = + tAcK^, 

respectively. Here, and tac are the amplitudes of the normal and Andreev reflections, 
respectively; tc and Iac are the amplitudes of the normal and Andreev transmissions, respec¬ 
tively. A, B, C, and D are the amplitudes of electronlike and holeslike quasiparticles in the 
SG region. All the amplitudes in Eq. ffTOj) can be determined by demanding wave function 
continuity at the interfaces. These boundary conditions are given by 


(0) = T 2 (0), T 2 (a) = T 3 (a). 


( 11 ) 


The scattering amplitudes can be obtained by numerically solving these continuity equations. 

The electron and hole operators of the outgoing states are related to the electron and 
hole operators of the incoming states via the scattering matrix—> 1 ^, 



where the element S'®® gives the outgoing electron current amplitude in response to an 
incoming electron current amplitude and S*®^ gives the outgoing hole current amplitude in 
response to an incoming electron current amplitude. The generalized current operator for 
the electrons and holes in the left electrode can be written as^^A^ 


h it) = 


27ih 




< at (E) au (E') >-<bl(E) {E') > 

+ <bt,(-E)b,A-E')> 


(13) 
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Substituting the scattering matrix in Eq. flT^ . we can obtain 

(E) A^P (a, E, E') {E') 

(i5)S;„„„(i5')07»p(B')J ’ 

where 

/l™” (a, E, E') = - Y, S2„„ (B) (E). 

n 

The general expression for the current fluctuations between contacts a and [3 is 
Sap = A4 (t) 5//3 {t') + 51/3 {t') A4 (t) >, 

with its Fourier transform 




p^rri'yp 


(14) 


(15) 


(16) 


27^6 {u — uj') Sap (w) =< A/q, (o;) 5Ip {uj') + 5Ip {uj') Ala (i^) > • (17) 


We restrict our consideration to coherent tunneling and neglect the Coulomb interaction. 
In the zero-frequency limit with a; = 0, providing all the information above, we can express 
the noise power as 


Sap = 



A^:;{a,E,E)A’^(^,E,E) 
x/„, (B) [1 - /,, (E)] + /.„ (E) [1 - (E)] 

+ E„, (E) (E) (E) (E) 

x/„, (-E) [1 - U (-E)] + /„ i-E) |1 - (-E)] 


By introducing the distribution function for the electrons fe {E) = 

[exp [{E — eV) /ksT] + 1]”^ and that for the holes fh (E) = [exp [{E + eV) /ksT] + 1]~\ 
the zero-temperature conductance can be obtained from the current operator and from the 
usual quantum statistical assumptions for the averages and correlations of the electron and 
hole operators in the normal reservoirs as 

G (eV) = Gq [ (^ — Wc'^ + Wac'^ -^ cos a da. (19) 

Jo V cos a / 

The shot noise power can be obtained as 


^ ( 0 ) = AeGo 


\rct \tc? + \rAct Jac 


cos a 


cos a 


cosada, 


( 20 ) 


where Go = 4e^A (eV) /h is the ballistic conductance of metallic graphene, V is the bias 
voltage, and N (e) = {e + Ep) w/ (Tihvp) denotes the number of available channels for a 
graphene sample of width w. 
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III. NUMERICAL RESULTS AND DISCUSSION 


Now we present numerical results for the Andreev reflection coefficients and the tunneling 
conductance for the G/SG/G junction with Uo 7 ^ 0. In this condition, there is a large 
mismatch of Fermi surfaces on the G and SG sides. Such a mismatch is well known to act 
as an effective barrier-'ii'I^. The transport properties of the two cases of Ep S> Aq {Ep = 0) 
and Ep <C Aq are significantly different. In the case of Ep 3> Aq, i.e., the incident electron 
and the reflected hole both lie in the conduction band, which results in the “retro-Andreev 
reflection”, while in the case of Ep -C Aq, only the “specular Andreev reflection” takes 
place, since the incident electron in the conduction band is converted into the reflected hole 
in the valence band^'^'^'^. In general, it is very difficult to reach the regime -C Aq in 
experiment^. So we only consider the case of S> Aq and the condition of comparable 
Ep and Aq, in the latter of which both normal Andreev reflection and specular Andreev 
reflection play roles and we can see that the retroreflection crosses over to specular Andreev 
reflection. 

Firstly, we consider dependence of the transport properties on the SG-barrier thickness. 
The tunneling conductance through the SG-barrier as a function of the thickness for different 
potential strengths Uq is shown in Fig. 2. Dimensionless thickness Fq® is used. The solid 
and dashed lines correspond to Uo/Ep = 2 and 10, respectively. When the bias voltage is 
small, the conductance exhibits oscillation features. Their oscillation amplitudes decay with 
increase of the thickness of the SG layer. The oscillation period of the shot noise power is the 
same as that of the conductance. The shot noise characterises correlations of the current. It 
can be seen from the curves that the growing trend of the shot noise is opposite to that of the 
conductance. When the conductance reaches the maximum value, the shot noise approaches 
the minimum value. This can be interpreted by the relation between the shot noise and the 
properties of the scatterer. A coherent conductor with all the transmission channels open 
(The open channel means that the transmission probability is close to one.) has minimal 
shot noise with the Fano factor approaching 0. A coherent conductor with all closed channels 
(The closed channel means that the transmission probability is close to zero.) has maximal 
shot noise with the Fano factor approaching 1. The strength of the shot noise is in the middle 
of the two limits when a conductor has open and closed channels coexistent. Therefore large 
transmission probabilities enhance the conductance and suppress the shot noise. It can also 
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be seen in the panels (b) and (d) that for the case of Ep = Aq the conductance and shot 
noise decrease with the increase of the thickness of the SG layer. This is because that the 
parameter k is proportional to (A/hni?) sin/? for identical Uq and it is in the exponential 
form of in the wave function as shown in Eq. ([8]). It is larger for the case of Ep = Aq 
than for the case of Ep S> Aq, which results in the quick decrease of the shot noise and 
the conductance. On the other hand, we can see from Fig. 2 that with the increase of the 
potential strength, the tunneling conductance decreased, which illustrates that the Fermi 
surface mismatch between the normal and superconducting regions suppresses transmission. 
We also considered the dependence of the tunneling conductance on the value of Uq for small 
bias voltages. As expected, we found that the oscillation amplitude decreases monotonically 
with the increase of Uq in the case oi Ep = Aq and hnally approaches a constant value. 

In Figs. 3 and 4, we provide numerical results of the tunneling conductance as a function 
of the bias voltage V. Similar results to Beenakker’si are obtained. In the limit of F —?■ 0, all 
the conductances have the same value of 4/3. For the case of Aq > Ep, a sharp change in the 
conductance occurs at eV = Aq and all the conductances vanish at eV = Ep, which is shown 
in Fig. 3. This is because of that no Andreev reflection occurs for all the incident angles 
(the critical angle of incidence ac = 0) when e = Ep. For the small bias voltages before 
the turning point eV = Ep, the conductance decreases with the increase of the SG energy 
gap. The conductance curves exhibit oscillatory behavior in the region of eV > Aq, which 
is different from the condition considered by Beenakker. By analyzing the transmission 
coefficients of the system, we found that in the condition of eV > Aq the transmission 
spectrum demonstrates oscillatory behavior and the oscillating period increases with the 
increase of the SG energy gap. The oscillatory behavior originates from the effect of the 
quantum-mechanical interference between the electron-like and hole-like quasiparticles^i^^^ 
in the SG barrier. This effect gives rise to oscillations in the reflection and transmission 
probabilities for the incident energies larger than the gap energy. 

In Fig. 4, numerical results of the conductance for different SG-barrier thickness a are 
provided. It can be seen that the conductance is sensitive to the SG-barrier thickness, 
especially in the condition of eV > Aq. In this condition, the oscillation period increases 
and the oscillation amplitude decreases sharply with the increase of the thickness. This is 
also a result of the exponential term in the wave function. It can be interpreted by 
the a —)■ cxo limit. In this limit, the components of the wave function with the exponential 


term are nonphysical, therefore only two of the components of wave function in the 
SG-region are physical. The model reduces to that of the the G/SG-junction, in other words, 
the model proposed by Beenakker- is obtained. 

We also considered the shot noise properties of the conductance through the G/SG/G 
structure. Numerical results of the shot noise and the Fano factor are provided in Fig. 
5. The parameters of panels (a) and (c) are the same with Fig. 4 and teh parameters of 
(b) and (b) are the same with Fig. 3. We can see from the curves that the shot noise 
oscillates greatly in the region of eV > Aq, while it increases monotonically in the region 
eV < Aq. These behaviors are similar to the tunneling conductance. In tunneling through 
the G/SG/G-structure, transmission is enhanced by the active hole channels. As a result, 
the values of the shot noise are small in comparison with the Poisson value, which is 2eG 
corresponding to the uncorrelated transport. When the hole channels in addition to the 
electron channels contribute to the transport, the interference effect is strong and the shot 
noise is signihcantly suppressed. In the region of eV < Aq, the impact of the proximity 
effect in the SG-region is strong giving rise to strong conductance and the values of S and F 
approach 0. In the case of eV > Aq, the amplitude of oscillation decreases with the increase 
of eV] in the case of eV < Aq, the amplitude of oscillation increases with the increase of 
eP, which originates from the same reason as the conductance. 


IV. CONCLUSIONS 

Based on the Dirac-Bogoliubov-de Gennes equation and the scattering theory, we inves¬ 
tigated the transport properties of the relativistic electrons and holes through the G/SG/G 
junction. We have deduced the analytical formulas of the tunneling conductance and the 
shot noise. Numerical results of the tunneling conductance and shot noise in the system 
are provided. We compared the two cases, one of which is in the presence of the specu¬ 
lar Andreev reflection and the other of which is in the absence of the specular Andreev 
reflection. The physical results can be summarized as follows. Firstly, the conductance 
increases with the increase of the thickness of the SG-layer and the shot noise is suppressed 
by the conductance in the case of eV = Ep- Secondly, the potential strength signihcantly 
affects the transport properties. It suppresses the conductance and enhances the shot noise. 
Thirdly, we obtained similar results with the model of the G/SG junction in the condition 
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eV < Aq. In the limit of a —)■ oo, the results of Beenakker for a G/SG junction can be 
reproduced. Fourthly, the thickness of the SG-layer affects the conductance and the shot 
noise more prominently in the condition eV > Aq, causing the decrease of the oscillation 
amplitudes and the characteristic features of the specular Andreev reflection. In conclusion, 
the conductance is a combined result of the Andreev reflection and the specular Andreev 
reflection, which can be tuned by the system parameters; the shot noise is suppressed by 
the SG-barrier because of the contribution of the hole channels in addition to the electron 
channels. 

V. ACKNOWLEDGEMENTS 

This project was supported by the National Natural Science Foundation of Ghina (No. 
11004063), and the Fundamental Research Funds for the Gentral Universities, SGUT (No. 
2014ZG0044). 


10 



^ C. W. J. Beenakker, Phys. Rev. Lett. 97, 067007 (2006). 

^ C. W. J. Beenakker, Rev. Mod. Phys. 80, 1337 (2008). 

^ A. H. Castro Neto, F. Guinea, N. M. R. Peres, K. S. Novoselov, and A. K. Geim, Rev. Mod. 
Phys. 81, 109 (2009). 

^ C. Bai, Y. Yang, X. Zhang, Phys. Rev. B 80, 235423 (2009). 

^ C. Bai, X. Zhang, Phys. Rev. B 76, 075430 (2007). 

® M. I. Katsnelson, K. S. Novoselov, and A. K. Geim, Nature 10, 1038 (2006). 

^ A. F. Andreev, Sov. Phys. JETP 19, 1228 (1964). 

® G. Bai, J. Wang, and Y. Yang, Superlattices and Microstructures 49, 151 (2011). 

^ G. Benjamin and J. K. Pachos, Phys. Rev. B 78, 235403 (2008). 

M. Titov and G. W. J. Beenakker, Phys. Rev. B 74, 041401(R) (2006). 

S. Bhattacharjee and K. Sengupta, Phys. Rev. Lett. 97, 217001 (2006). 

G. Bai, Y. Yang, Phys. Lett. A 374, 882 (2010). 

M. J. M. de Jong, G. W. J. Beenakker, Phys. Rev. B 49, 22 (1994). 

G. W. J. Beenakker, Phys. Rev. B 46, 19 (1992). 

G. E. Blonder, M. Tinkham, and T. M. Klapwijk, Phys. Rev. B 25, 4515 (1982). 

Ya. M. Blanter, M. Biittiker, Phys. Rep. 336, 1 (2000). 

B. A. Muzykantskii and D. E. Khmelnitskii, Phys. Rev. B 50, 6 (1994). 

M. P. Anantram, S. Datta, Phys. Rev. B 53, 16390 (1996). 

G. X. Bai, X. D. Zhang, Phys. Lett. A 372, 725 (2008). 

G. X. Bai, Y. L. Yang, Appl. Phys. Lett. 92, 102513 (2008). 

W. L. McMillanand, P. W. Anderson, Phys. Rev. Lett 16, 85 (1966). 

K. Pasanai, Physica G 499, 1 (2014). 

Z. C. Dong, R. Shen, Z. M. Zheng, and D. Y. Xing, Z. D. Wang, Phys. Rev. B 67, 134515 
(2003). 


11 




FIG. 1: Schematic illustration of reflection and transmission of quasipaticles in the 

graphene/super conductive graphene/graphene (G/SG/G) heterostructure, which occupies the x-y 
plane. The left G region extends from x = —oo to x = 0 and the right G region extends from x = a 
to X = + 00 . The superconductive region occupies 0 < x < a. A voltage is applied between the left 
and right G regions. As a result of scattering by the SG, an incident electron can be transmitted 
and reflected both as an electron and as a hole into the right and left electrodes, respectively. The 
bottom is the superconducting potential profile with a the SG-region width and A the value of the 
superconducting gap. 
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FIG. 2: Tunneling conductance and shot noise of the G/SG/G junction as a function of the SG- 
thickness k^a for different potential strengths. Solid and dashed lines correspond to the potential 
strengths of Uo/Ep = 2 and 10, respectively. The other parameters are 0 = 0, s/Ep = 0.005, and 
a = 5. 


13 

























eV/A 


FIG. 3; Tunneling conductance as a function of the bias voltage for different energy gaps of 
the SG. The solid, dashed, and dotted lines correspond to the energy gap in the SG region of 
Aq/Ef = 1, 2, and 5, respectively. The other parameters are </> = 0, Uq/Ef = 100, and a = 5. 
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FIG. 4: Tunneling conductance as a function of the bias voltage for different thicknesses of the SG 
region. The solid, dashed, and dotted lines correspond to the SG-region thickness a = 20, 30, and 
50, respectively. The other parameters are </> = 0, U^/Ep = 10, and Aq/Ep = 0.1. 
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FIG. 5: Shot noise and the Fano factor as a function of the bias voltage for different energy gaps 


of the SG. The parameters in panels (a) and (c) are the same as Fig. 4. The parameters in panels 
(b) and (d) are the same as Fig. 3. 
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